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We perform a non-perturbative study of the Coleman- Weinberg phase 
transition in scalar QED. Our method permits a consistent treatment of the 



5-H ' effective potential near the origin, a region not accessible to perturbation the- 

ory. As a result, we establish reliably the first order character of the phase 
transition for arbitrary values of the scalar coupling. 
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The seminal Coleman- Weinberg one-loop computation ||T| of the effective poten- 
tial in scalar quantum electrodynamics revealed the phenomenon of spontaneous 
symmetry breaking induced by quantum fluctuations. This perturbative calculation 
indicates a first order transition between the symmetric phase and the one with 
spontaneous symmetry breaking in dependence on the scalar mass term. This ob- 
servation has motivated many theoretical investigations, ranging from a lower bound 
on the Higgs mass in the standard model (for top quark mass below 80 GeV) - the 
Linde- Weinberg bound [Q] - to speculations about the first order character of the 
electroweak phase transition at high temperature 0. The latter could have impor- 
tant consequences for the creation of the baryon asymmetry in the early universe 
01 by providing a new scenario for non-equilibrium conditions at temperatures of 
the order of the Fermi scale. Another area of applications are models of inflationary 
cosmology where Coleman- Weinberg type potentials are frequently used 0. 

Strictly speaking, however, the perturbative calculation is only valid for the 
convex part of the effective scalar potential where all mass terms for fluctuations 
are positive or zero and the saddlepoint expansion is well defined. The failure of the 
loop expansion for negative mass terms for the fluctuations manifests itself through 
an imaginary part in the effective potential. The region in field space relevant for the 
first order character of the phase transition does not belong to the convex part of the 
potential and is therefore not directly accessible to perturbation theory. Since the 
problems are related to the fluctuations of scalar fields and the scalar interactions 
are small in the critical region, one usually discards the scalar fluctuations and only 
considers contributions from gauge fields. Indeed, the quartic scalar coupling A is 
much smaller than e^ (with e the abelian gauge coupling) in the critical region and 
the neglection of scalars seems justified at first sight. Unfortunately, this assumption 
cannot be verified within the standard perturbative framework since a saddle point 
expansion is crucial for its validity. Furthermore, little is known about the effect of 
higher-order loop corrections on the weak first order character of the transition. 

Recently, a new non-perturbative method for the study of spontaneous symmetry 
breaking has been proposed. It is based on the concept of the average action ^ 
which is the continuum analogue of the blockspin action [0 in lattice theories. The 
average action F^ is the effective action for averages of fields taken over a (Euclidean) 
volume ~ k~^. It is obtained by integrating out the quantum fluctuations with 



momenta q^ > k"^. The dependence of F^ on the "average scale" k is governed by 
an exact non-perturbative evolution equation [§] 0. It can be shown f^ that the 
effective average action interpolates between the classical action as defined at some 
short distance scale A~^ (Fa = S) and the generating functional of the IPI Green 
functions - the usual effective action F - for k = (Fq = F). The properties of the 
phase transition can be read off from the average potential Uk which is a manifestly 
real quantity for all values of the scalar field. Even for regions where perturbation 
theory ceases to be valid the non-perturbative evolution equation - which takes the 
form of a renormalization group- improved one-loop flow equation - still applies. This 
method seems therefore suitable for an investigation of non-perturbative aspects of 
the Coleman- Weinberg phase transition. 

The concept of the average action was generalized for abelian gauge theories ||10[ , 



nH . The dependence of the effective scalar potential Uk [|TOl and the gauge coupling 



Tl[| on the average scale k was computed in arbitrary dimensions, demonstrating 
how the new non-perturbative method can successfully cope with the infrared diver- 
gences present in the standard perturbative approach. The exact non-perturbative 
evolution equation for the scale dependence of the effective average action has been 
adapted to non-abelian [|12] and abelian |jl3[ gauge theorie^. It appeared that the 



flow equations for the scalar potential and the gauge coupling proposed in |]IO[,||TT 
coincide with the exact non-perturbative evolution equation in the limit where the 
effective action is truncated to include only standard kinetic terms (with wave func- 
tion renormalization) for the scalar and gauge fields as well as a scalar potential. 

In this letter we employ the evolution equations for a non-perturbative study of 
the phase transition of the abelian Higgs model. We restrict the discussion to four 
dimensions and vanishing temperature and study the phase transition as a function 
of the mass term, as in the original papers fl]], |0|. We are mainly interested in 
the phase structure of the theory and concentrate on the scale- dependent effective 
scalar potential Uk{p), where p = (p*(p. The exact evolution equations for Uk can 
be infered from ref. [|l3l, but are far too complicated to be solved exactly. We 
therefore proceed to different approximations - all of them non-perturbative and 
going substantially beyond the one-loop computation. We first consider the limit 



'''See ref. |^ for earlier versions of exact renormalization group equations. 
^See ref. MM for a different approach. 



where the scalar fluctuations can be neglected and deal with the scalar fluctuations 
later. Furthermore, we approximate the derivative terms in F^ by standard kinetic 
terms for both the scalar field and the gauge field and neglect the wave function 
renormalization of the scalar field. The flow equation for the scale dependence of 
the average potential becomes very simple in this case and reads, with t = \n{k/A), 

dtUM = --^e^ / dx xd,{P{x) + 2e^p)-\ (1) 

OTT^ JO 

Here primes denote derivatives with respect to p, e^{k) is the running abelian charge, 
and the effective inverse propagator 

P{x) = x l-exp('-^') (2) 

is a function of the momentum squared q^q^ = x, which contains an effective infrared 
cutoff -P(O) = k^. The partial derivative dt on the r.h.s. of eq. (|1|) is meant to act 
only on P. 

After performing the momentum integration the flow equation (|l]) becomes a 
partial differential equation for Uk as a function of the two variables p and k. It has 
to be supplemented by an equation for the running of the abelian gauge coupling 
e^{k). Since the latter is only logarithmic - the usual one- loop (3 function is a good 
approximation (cf. eq. ([T0|) ) - we may first simplify for constant e^. The evolution 
equation ([l|) can then be solved easily in the form 

Uk{p) = ^,e'J^ dxx{P + 2e'p)-\ (3) 

Taking the limit fc — > 0, we obtain the effective potential (apart from an irrelevant 
constant) 

?, . . r / n \ 1 1 

(4) 



f/o(p) = m^p + -^^e p 



,n U-\ - I 
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where po is fixed through Uq{pq) = 0. This is exactly the well-known one-loop result 
of Coleman and Weinberg 0. The mass term at the origin ml is obtained as an 
integration constant and may hence be treated as a free parameter. The second 
free integration constant corresponds to the quartic scalar coupling taken at some 
appropriate scale, i.e. U"{p = p^), or, equivalently, the scale po. As a function of 
ml the theory indeed exhibits a first order phase transition in this approximation! 
Beyond the approximation of constant e^ the logarithmic running of e^(/c) does 
not alter the qualitative features of the phase transition. A numerical solution of 



eq. (0) with running e^(/c) gives only minor quantitative corrections. This also holds 
if the gauge boson contribution to the wave function renormalization (anomalous 
dimension) of the scalar field is included (see below). The only remaining difference 
between the exact evolution equation and the flow equation in the limit of neglected 
scalar fluctuations concerns then a generalization of the gauge boson kinetic term in 
Ffc beyond the standard form ~ F^^F^^ and a similar generalization of the covariant 
scalar kinetic term. These effects would modify the truncated form (-P(x) + 2e^p)~^ 
of the effective gauge boson propagator in the presence of a constant background 
scalar field p. In principle, the scale dependence of these modifications can in turn 
be determined from exact flow equations. In the limit of small gauge couplings they 
could alternatively be approximated by one-loop expressions - this would allow a 
determination of the e^-terms in the /5-functions without ever doing a two-loop cal- 
culation, in complete analogy to the successful demonstration in pure scalar theories 



| T5[| . It is easy to convince oneself that none of these small propagator modifica- 
tions could possibly alter the first order character of the phase transition indicated 
by the flow equation (1). The dominant effect — U'^{p) turning negative for small 
enough p as k decreases - is independent of the details of the effective propagator. 
We conclude that in the limit of neglected scalar fluctuations the form of the av- 
erage potential Uk indeed indicates a first order phase transition. For small e^ the 
Coleman-Weinberg potential (^ gives an accurate quantitative description, with 
only minor modifications ~ e^ as compared to a more complete non-perturbative 
analysis. 

The crucial question in order to firmly establish the first order phase transition 
concerns, however, the contributions to the effective potential from fluctuations of 
the scalar field. Are they small as suggested by the small ratio A/e^? In order to 
give a non-perturbative answer to this question we truncate the average action (with 
background gauge fixing [|1^) 

In this approximation and for a — i> the flow equation for the average potential 
reads 



("OO 

dtUk{p) = 3v4 dx xdtln[P + 2e'^p] 
Jo 



+Vi / dx xdt \n[P + U'^ip) + 2pU'^{p)] 
Jo 

+vj dx xdt\n[P + U'kip)], (6) 

Jo 

with Vi = l/327r^. As in (1), the partial derivative dt acts only on P. We have 
solved this equation numerically and the corresponding phase diagram is depicted 
in fig. 1. The basic variables needed in order to describe the shape of the aver- 
age potential are the location of the minimum po{k) (where f/((po(^)) = 0) and 
the quartic coupling related to f/^'(po)- It is convenient to use dimensionless renor- 
malized parameters, in particular e^(A;), the renormalized abelian charge, A(A;), the 
renormalized quartic scalar self-coupling, and K,{k), the dimensionless renormalized 
location of the potential minimum. Their definitions are as follows: 

e'{k) = Z,lf{k) (7) 

K{k) = k'^Z^^kPoik) (8) 

A(A;) = Z-p'^Po), (9) 

where Zip^k{ZF,k) denotes the wavefunction renormalization of the scalar (gauge) 
field. In the projection of the phase diagram on the A, k, plane we observe a crit- 
ical line separating the symmetric phase (SYM) from the phase with spontaneous 
symmetry breaking (SSB). The other curves represent the fiow of K,{k) and A(A;) 
for initial values /t(A) and A(A) given at a fixed short distance scale A. (We use 
e^(A) = 0.1 .) The arrows represent the direction of the fiow with A; ^ 0. Tra- 
jectories starting below the critical line end with k, = and, therefore, with the 
potential minimum at the origin. For trajectories above the critical line k, finally 
diverges ~ k~'^ such that po reaches a constant value corresponding to the square of 
the vacuum expectation value of ip. 

For A > 10~^ all quantitative details of the numerical solution can be understood 
by a simple polynomial expansion of Uk around po{k), setting U'j!'{po) = and 
similarly for higher derivatives. In this approximation the evolution equations for 
the SSB regime read 

9te' = ^,e's',{2X^,2e'^) (10) 

dtX = -^^\\9st{2\K) + l) + ^e'st{2e\) 

-—e^Xsl^{2e^K, 2Xk) + —X^Kml^{2XK, 0) (11) 



+ ^e^/ts^ j2e^K, 2Ak) - -^X'^K'^mlJ2XK, 0). (12) 

(The last two terms in eqs. (pH]) and ([T^) arise from the anomalous dimension 
rjip = —dt In Zfc,<^ and correspond to 2?7^A and —rj^pK, respectively.) The new elements, 
compared to the standard perturbative one-loop /3-functions, are the threshold func- 
tions s^, sf, s|, s-j^^ and 1712^2- They describe the effective decoupling of particles with 
mass much larger than k. Indeed, two different masses play a role in the regime 
with spontaneous symmetry breaking (SSB regime). The mass of the gauge boson 
at the scale k is given by 

M\k) = 2e\k)po{k)Z^Zp' 

= 2e\k)K{k)k^ (13) 

whereas the mass of the Higgs boson corresponds to 

m\k) = 2po{k)U'j:{po)Z-' 

= 2X{k)K{k)P. (14) 

The physical masses of these particles (as defined by the two-point function at zero 
momentum) are obtained from eqs. ([l3|) , ([l4|) in the limit k ^ 0. The threshold 
functions are all normalized to one in the limit of vanishing arguments. For k, — > 
the evolution equations (|10D,(|lT]) therefore reduce to the standard perturbative 
renormalization group equations. On the other hand, the threshold functions vanish 
for large arguments M^/fc^ or rn^ /k'^. Their precise form depends on the precise 
form of the infrared cutoff function Rk{x) (for details see ref. [|l3])- With Rk{x) = 
Z^_fca;exp(— p-)[l — exp(— p-)]^-*^ they read for large values of the argument {w ^ 1) 

[17]4(w;) = 6nC(3)ur("+i) 

-12n{n + 1) [C(3) + C(4)]it;"("+') + 0{w-^''+^^) (15) 

and similarly for the other threshold functions. The appearance of threshold func- 
tions in the /5-functions is physically expected. Particles with masses much larger 
than k should not influence the running of couplings with the infrared cutoff scale k. 
Their mass already acts as an effective infrared cutoff and the variation of a second, 
much smaller cutoff k should only have a small effect. 



The second new ingredient is the evolution equation for k. It is closely related 
to the "quadratic renormalization" ^j oi po, or, equivalently, the scalar mass term. 
It describes within a renormalization group framework the physics related to the 
"quadratic divergences" appearing in perturbation theory. Such an equation is nec- 
essary in order to give a meaningful description of the threshold effects since those 
must depend on the ratio of the renormalized (fc-dependent) mass over k. Although 
well known in practice in the Wilson approach to renormalization (for example 
lattice studies), the effect of quadratic renormalization is missing in many versions 
of the renormalization group equations. This happens typically if there is no in- 
dependent running scale to which the masses can be compared, as for example in 
the Coleman- Weinberg analysis of the effective potential were both the gauge boson 
mass and the running scale are given by the "classical" or "background" field cp. 

In our approach the renormalization group equation for the fc-dependent mini- 
mum of the average potential arises naturally from the minimum condition 

JtU'kiPoik)) = ^fiiPo) + Ul{pof-^ = (16) 

and the evolution equation for i9(f/^ which involves terms ~ e^k"^ and Afc^. This also 
explains the perhaps somewhat surprising negative power of A in the third term of 



eq. (12): For very small A even a moderate change in f/^((9tf/{, ~ e^fc^) induces a 
large change in the location of the minimum. 

Given a set of initial values e^(A), fi;(A), A (A) at some short- distance scale k = A, 
the evolution equations describe how the various couplings run in dependence on k. 
The solution for A; — > gives the vacuum expectation value pQ{k = 0) as well as the 
renormalized gauge coupling and quartic coupling at zero momentum. The value 
of po{k = 0) determines whether the theory is in the symmetric phase (po(0) = 0) 
or in the SSB phase (po(0) > 0). For a theory in the symmetric phase k reaches 
zero at some scale ks > 0. Once K,{ks) = 0, the evolution has to be continued 
using equations appropriate for the symmetric regime. Since most of the physical 
properties can be derived from the equations in the SSB regime, we will not present 
these equations explicitly here and refer the reader to ref. [|l^. For a theory in 
the SSB phase k remains positive for all values of k and eqs. (p!0D-(p!2D are sufficient. 
This system of coupled non-linear differential equations still looks quite complicated. 
It allows, however, for analytical solutions in some approximations which make the 
main characteristics easy to understand. 



First we discuss the linear regime defined by 2Ak ^ 1, 2e^fi; <^ 1. All masses are 
then much smaller than the scale k and we should recover the results of perturbation 
theory. The evolution equations simplify considerably: 

dte' = ^v,e' (17) 



dA = 24f4 



o 



dtK = -2k + 8f4 - I2V4XK + 12t;4— - 24t;4-— + 12f4e\. (19) 

A A 

Equations (0), ( |l8|) coincide with the standard one-loop result in perturbation 
theory. The third and fifth term of eq. ([19D are obtained through the expansion 
of sf{w) for small w: sf{w) = 1 — w + 0{w'^). Combining eqs. (|18D and ( p^QJ) one 
can obtain the anomalous mass dimension of the scalar field in one-loop order. The 
running of the gauge coupling 

e\k) = ,^/^^^,,, A (20) 

1 + 3^746^ (A) In -^ 

is a small logarithmic effect and will be neglected in the following (e^(A;) = e^). 
There is no fixpoint for the running of A (or for the evolution of the ratio A/e^) 
and the positive (3 function drives A always towards zero. The interesting physics 
which is associated with the first order character of the Coleman- Weinberg-type 
phase transition occurs for values of the quartic coupling A much smaller than e^, 
typically A ~ v^e"^. For A <^ e^ we can neglect the terms ~ A^ and e^A in eq. (|1^). 
The running of A is then easily found 

Xik) = 2Av4eHn(y\ (21) 



fcc = Aexp(-,^^ . (22) 



A(A) 

24^46^ 

According to eq. (|18[) A reaches zero at some scale k^ ^ and becomes negative for 
k < kc- To solve eq. (^) for k, in the limit |A| ^ e^ we observe the existence of an 
UV-stable fixpoint for Xk: 

dtiXK) = -2[XK-{XK),] (23) 

(Xk), = Gv^e^ (24) 

8 



This fixpoint is represented by the dashed hne in the phase diagram in fig. 1. Com- 



bining the solution of eq. (^) with eq. (^) we find 



A2 

I I — 1~ / \ A 



1 I + Aa^ 



4e2 \n{k/k,) 

We have parametrized the initial value at A such that Aa measures the deviation 
of Xk, from its fixpoint value of eq. (^) at the scale A. 

Due to the running of A there is no fixpoint in k. For Aa = one sees that n 
diverges for k -^ k^. On the other hand, for sufficiently negative Aa, K,{k) will vanish 
at some scale kg > kc and the evolution will continue in the symmetric regime. For 
k < ks we replace the parameter K,{k) by 

m'M = Z-p'M (26) 

and use evolution equations adapted to the symmetric regime. In particular, the 
scale dependence of the mass term simplifies for |A| ^ e^ to 

dml = -^,e'k\ (27) 

The mass term at A; = is therefore connected to kg by 

"• = i^.''^i- (28) 

There is a critical value 

separating these two different behaviours. For Aa > Aa^ we see that K{k — ^ kc) 
diverges, whereas for Aa < Aa^ we end in the symmetric phase. For Aa = Aa^ we 
find that Hi{k —>■ kc) reaches a finite non-zero value 

1 
2e 

Thus the system of equations ([l7[) -([l9[) exhibits a bifurcation phenomenon at the 
critical value Aac: No finite non-vanishing value of K,{kc) except 1/26^ can be 
reached. 



(30) 



In the region of initial values Aa < SAa^ the linear approximation of eqs. (|T 
j) remains always valid. In this region the mass term is bounded from below 



,2^ 3 0,2 _ 3 ^2a2„__/ -^(^) 



In the SSB phase k always diverges and there must be some regime where the linear 
approximation breaks down. A quick inspection of the full equations ([TTl) (|12D shows, 
however, that for A > -^e.^ the threshold functions only imply a stop of the running 
of A and po for 2e^fi; > 1. We can therefore approximate 

9iA = 1 f 2e'^K > 1 

for S 32 

5tPo = J I A > sf,e^ 

For A A > we find that 2e^K becomes larger than one at a scale kp where 
\n{k'^p/k1) > 1, A(A;i?) > of^e^- For this region there is a lower bound on the vacuum 
expectation value and the Higgs mass 



2e2 '^ 2e2 ^ \ 2Avie* J 

m > — — — — e k^. (33) 



2 ^ 3exp(l)^2,.2 

r2 



i-K^ 



In order to make a quantitative comparison with the perturbative Coleman- 
Weinberg analysis |1| for A = O(e^) we want to identify the value /t(A) which leads 
to the "massless theory", i.e. m1{k = 0) = 0. For a quartic potential at the short 
distance scale A (and Z^^a = 1) we observe 

ml(A) = -A(A)fi:(A)A2. (34) 

Using the flow equations (^) and (^) allows us to express the initial deviation Aa 
from the (Afi;)-flxpoint (24) in terms of the mass term at the origin m^ for k = 0: 

Therefore, the "massless theory" corresponds to Aa = 0. Using the linear equations 
(0),(^) as long as 2e^fi; < 1 and eq. (|32D for 2e^fi; > 1 we deduce for m^ = the 
mass relation (for A; = 0) 

(36) 



m2 _ A(A;;,) _ 3 ^ 

'M^ ~ c2 ~ R^*^ ' 



e" avr" 



10 



where kb is determined by 2e^K.{kh) = 1. Thus, in this approximation the well-known 
result of Coleman and Weinberg |Q is recovered. 

We conclude that except for the immediate vicinity of Aa^(2Aa^ < Aa < 0) 
the system of equations (|10D-(|12]) reproduces the results of the Coleman- Weinberg 
analysis, including the effect of dimensional transmutation (cf. eqs. (|3TD,(|33|) and 
(il))- I^ is tempting to think that the bifurcation phenomenon observed in the 
linear approximation corresponds to a first order transition, leading to mass bounds 
somewhat weaker than (^) and (^). We observe that the bifurcation occurs for 
values A <^ e^. If the effect of scalar fluctuations could be neglected for very small 
A, we may use eq. (|l|) to establish the first order character of the phase transition. 
Using the system of flow equations ([I0|)-(|r2|) for A > e^ and eq. (|l|) for A < e'' would 
then give a detailed non-perturbative quantitative picture of the Coleman- Weinberg 
phase transition]^. 

For a numerical study of the evolution equation (^ for A ^ e^ we have to con- 
sider potentials which may have more than one local minimum. For this reason we 
use for our numerical investigation an approximation of Uk as a quartic polynomial 
in p, i.e. we include couplings up to ip^. As long as Uk exhibits a local "asymmetric 
minimum" (A > 0, k > 0), the effective average potential is parametrized in terms 
of the four variables A, K,ml/k'^, and A^, where Xsik) = Z^']jJ',^{p = 0) denotes the 
renormalized quartic scalar self-coupling at the origin. This parametrization is best 
suited for the study of a (weak) first order phase transition: We expect that the sys- 
tem near the phase transition develops a small (positive) mass term at the origin. 
Therefore, it is most important that all the relevant quantities like masses and quar- 
tic couplings are followed by their appropriate evolution equation. Higher couplings 
like f/("(po) appearing on the r.h.s. of the evolution equation for k and A are deter- 
mined as functions of the four variables using the polynomial approximation for Uk- 
Our parametrization yields information on the global structure of the effective po- 
tential since it contains information on the symmetric extremum (m^, A^) as well as 
on the asymmetric one (k. A). When the broken phase disappears, (A < or k = 0), 
the evolution is continued in the symmetric regime with the corresponding evolution 
equations. A more detailed account on this method and some generalisations will 
be given elsewhere [|T9|. 



^The two systems of flow equations have a smooth overlap for small e^. 
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The phase diagram in the range of small values of X{k) is shown in fig. 2. We 
have specified the initial conditions at some fixed short distance scale A by the 
requirement that the classical potential U/^ is quartic in the fields, U{p) = |A(A)(p — 
Po(A))^. For all trajectories, we fix e^(A) and A(A). The phase diagram is then 
obtained through the variation of k(A) which serves as a free parameter. 

The phase transition is clearly of the weak first order type and we observe three 
characteristic regions in the (A, k) plane: In region I the average potential has only 
one minimum. (This means m^ < for k, > 0.) For trajectories remaining within 
region I this holds for all values of k. Once K{ks) becomes zero, region I should be 
continued in the symmetric regime with the only minimum now at the origin. The 
behaviour of trajectories remaining within region I is very similar to a second order 
phase transition. In particular this concerns the values of A shown in fig. 1. For 
not too small values of A (A) the phase transition is therefore "almost of the second 
order" . The weak first order character of the transition is related to trajectories 
so close to the critical line that they cross the dashed boundary between regions 
I and II. For regions II and III the average potential has two local minima. The 
"asymmetric minimum" occurs for k > and is the absolute minimum in region 
II. The "symmetric minimum" is at the origin. This becomes the lowest minimum 
in region III. The (dashed) borderline separating I and II corresponds to mf = 0. 
Therefore the "massless theory" (the Coleman- Weinberg case) is described by the 
trajectory A for which m1{k = 0) = 0. On this trajectory the mass relation ( ^6|) 
is found to coincide with the numerical result within 4%. The regions II and III 
are separated by the dashed-dotted line. This borderline describes the effective 
potential with degenerate minima f/fc(0) = Uk{fo)- The trajectory with initial value 
Kcrit is running asymptotically towards a potential with degenerate minima and 
corresponds to the critical line of the phase diagram. All trajectories to the right 
of B with k{A) > Kcrit belong to the phase with spontaneous symmetry breaking, 
whereas all trajectories to the left with k{A) < Hcrit will end up for A; ^ with lowest 
minimum at the origin and therefore belong to the symmetric phase. The trajectory 
C is the spinodal line of the model. It describes, similarly to A, an asymptotically 
massless theory but now for scalar fields at the asymmetric local minimum. The 
first order part of the phase diagram can only be seen in the small interval of initial 
values ka with k^ < k/^ < k^. This corresponds to A smaller than approximately 
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22^e^ in the final part of tlie running. 

In conclusion, the weak first order character of the phase transition has now been 
established on firm non-perturbative grounds. In particular, we emphasize that our 
non-perturbative evolution equation (P) is valid for all values of p, in contrast to 
the perturbative treatment which becomes unreliable around the origin, as pointed 
out in the original work [^. The average potential is always real. The perturbative 
imaginary part has to be considered as a pure artefact of perturbation theory. The 
average potential becomes convex for k -^ as can be seen on very general grounds 
0,[H. This approach to convexity is described by the exact evolution equation, but 
any polynomial approximation for the potential breaks down in the SSB phase[[ We 
are not interested here in the physics related to the approach to convexity and stop 
the running of Uk once k^ becomes of the order of the curvature of the potential 
at a local maximum. At these values of k the masses m1{k),m^{k) or M'^{k) have 
already reached almost constant values (and similarly for the interactions), and we 
use in practice these values instead of the values for fc = 0. We can then proceed 
to compute the physical photon mass M^ (which is proportional to the vacuum 
expectation value po), the scalar mass in the asymmetric minimum m^, and at 
the origin ml as a function of the initial value k{A). As discussed before we may 
use ml to label the distance from the phase transition. In fig. 3 we have plotted 
the various masses as a function of m^ and indicated values corresponding to the 
curves A, B and C in fig. 2. The mass ratio m? /M"^ is given in fig. 4 together with 
the Coleman- Weinberg prediction (|36|) for ml = 0. The visible small difference 
between our curve and the prediction (^) is actually not a real effect due to the 
scalar fluctuations neglected for eq. (|36|) but rather an artefact of the polynomial 
approximation to the potential. This can be seen by using the numerical solution 
with polynomial approximation, where the scalar fluctuations are discarded. The 
resulting curve for m? /M"^ is essentially identical to fig. 4. We conclude that the 
effects of scalar fluctuations are really tiny for the small values of A relevant for the 
first order character of the transition. Our non-perturbative approach permits a very 
precise quantiative understanding of all physically interesting masses and couplings 
as well as of the complete shape of the effective potential. This is achieved through 
the combination of the numerical solution of eq. (^ in the range of A(A;) > lOe^ 

''For a discussion of the approach to convexity for the pure scalar theory see ref. Q . 
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(cf. fig.l) with the analytical solution (^ for A(A;) < lOe^. 

In summary, we find a first order phase transition for scalar QED for arbitrary 
values of the quartic scalar coupling A. The mass gap between the symmetric phase 
and the one with spontaneous symmetry breaking is exponentially small (except 
for tiny values of A). This is an aspect of dimensional transmutation ||l|] due to 
the logarithmic running of A. In the vicinity of the phase transition the running 
quartic scalar coupling is always very small. As a result, we find that the one-loop 
calculation of the effective potential in the abelian Higgs model |jl| gives a rather 
accurate description of all physical quantities of interest in the vicinity of the phase 
transition. Our non-perturbative method permits (for the first time) a consistent 
treatment of the region around cp = in the phase with spontaneous symmetry 
breaking. Comparison with our non-perturbative results gives a simple prescription 
of how to deal, in this particular situation, with the effects of scalar fiuctuations in 
the one-loop calculation: The unphysical imaginary part should be discarded or the 
scalar fiuctuations should simply be omitted, in the range of (f for which a negative 
mass squared for the scalar fiuctuations is predicted by the one-loop calculation. 
The one-loop potential computed with this prescription is very close to the average 
potential Uk at a scale kcv which is given by k"^^ ~ 2m^„^., with m^„^. = —U\pmax) 
corresponding to the curvature at the local maximum of the potential. 

More generally, the effective potential Uq is not the most appropriate quantity 
for a description of spontaneous symmetry breaking. On general grounds Uq is con- 
vex and therefore more subtle to handle in the phase with spontaneous symmetry 
breaking. On the other hand, the average potential Uk is not convex and the quan- 
tum effects leading to convexity of the effective potential Uq only set in for k < kcv 
TsU - As long as the relevant physical distances are smaller than k~^ the fluctuations 



leading to convexity can be neglected. This also holds for all properties related to 
the true vacuum (corresponding to the absolute minimum of the potential), such as 
the masses of physical particles and their interactions. We should mention, however, 
that the flattening of the potential due to the approach to convexity may play a role 
in certain cosmological scenarios of inflation. 

The present work has also permitted precise quantitative tests for analytical 
and numerical approximate solutions of the non-perturbative evolution equations. 
These evolution equations remain valid in three dimensions or for a non-vanishing 
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temperature where perturbation theory fails due to severe infrared problems. In 
the presence of large scalar couplings the one-loop approximation neglecting scalar 
fields is not expected to give a qualitatively correct description. The present work is 
therefore an excellent starting point for a non-perturbative treatment of the three- 
dimensional abelian Higgs model - which should describe the phase transition of 
superconductors - or of the high temperature phase transition in four-dimensional 
scalar quantum electrodynamics. 
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Figure Captions 

Fig. 1: The phase diagram of scalar QED. The thick hne separates the symmetric 
phase (SYM) from the phase with spontaneous symmetry breaking (SSB). 

Fig. 2: The phase diagram for small A. We indicate the "massless" (A), critical 
(B), and spinodal (C) trajectory. 

Fig. 3: The mass of the gauge field M as a function of the scalar mass term at 
the origin m^. We also indicate the scalar mass m for spontaneous symmetry 
breaking. All masses are given in units of the cutoff A. 

Fig. 4: The ratio of the Higgs-boson to gauge-boson mass as a function of m^. The 
star denotes the Coleman- Weinberg result. 
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This figure "figl-l.png" is available in "png" format from: 
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